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Criticality of an isotropic-to-smectic transition 
induced by anisotropic quenched disorder 
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We report combined optical birefringence and neutron scattering measurements on the liquid 
crystal 12CB nanoconfined in mesoporous silicon layers. This liquid crystal exhibits strong nematic- 
smectic coupling responsible for a discontinuous isotropic-to-smectic phase transition in the bulk 
state. Confined in porous silicon, 12CB is subjected to strong anisotropic quenched disorder; a short- 
ranged smectic state evolves out of a paranematic phase. This transformation appears continuous, 
losing its bulk first order character. This contrasts with previously reported observations on liquid 
crystals under isotropic quenched disorder. In the low temperature phase, both orientational and 
translational order parameters obey the same power-law. 

PACS numbers: 64.70.Nd, 61.30. Eb 



Since its introduction in 1975 by Imry and Ma [T], 
the influence of random fields on phase transitions has 
been one of the most debated topics in condensed mat- 
ter physics. Various realizations of this phenomenology 
could be achieved by confining liquid crystals (LCs) in ge- 
ometrically disordered porous materials or gels [2,-i4!. In 
principle, the geometric restriction introduces two forms 
of disorder in LCs: random orientational fields that cou- 
ple to the nematic director and random positional fields 
that couple to the smectic order. Studies on the influence 
of these couplings while varying the magnitude of disor- 
der and its isotropic or anisotropic character has largely 
contributed to recent progress in the understanding of 
quenched disorder (QD) effects on phase transitions in 
LCs m HJ [B]. Owing to their generic character, these 
flndings are also relevant to many other systems, most 
prominently for supraconductors or for superfluids [7]. 

Both the second order nematic-smectic A (N-SmA) 
and normal-superconducting transitions can be mapped 
onto each other and, in principle, fall in the universal- 
ity class of the 3D XY model [71 [8] , since both can be 
described by a complex order parameter representing the 
amplitude and phase of a sinusoidal- varying smectic mass 
density wave or a macroscopic wave function, respec- 
tively. However, as was pointed out flrst by de Gennes 
the coupling between nematic (Q) and smectic (77) order 
parameters (OPs) actually takes the N-SmA transition 
away from that class. Conversely, increasing the strength 
of isotropic QD can gradually shift the character of the N- 
SmA transition from tricritical back to 3D XY universal- 
ity. This was verified by studies on n-octyl-cyanobiphenyl 
(8CB) in aerogels or loaded with random dispersions of 
aerosil particles [3l H] . For these systems, the influence 
of strong isotropic QD in the weak Q-r] coupling limit 



was addressed. For anisotropic QD, this striking effect is 
already visible in the limit of weak disorder strength in 
the case of 8CB. Following the argumentation of Garland 
and lannachione [1], this suggests that the Q-rj coupling 
is reduced by isotropic QD, and is even entirely turned 
off by anisotropic QD. 

An interesting test of this hypothesis was made by Ra- 
mazoglu et. al. [8 on the lOCB/aerosil composite. In 
the bulk, this LC shows a strong Q-rj coupling, which 
leads to a direct and discontinuous isotropic (I) to SmA 
transition. Under isotropic QD, the transition remains 
first order but the temperature scaling associated to the 
growth of the order parameter tends towards tricritical 
behavior as the density of disorder is increased [8., 9, . The 
emerging question is whether a direct I-SmA transition 
of a strongly Q-rj coupled LC could become continuous 
under strong anisotropic quenched disorder. 

Here we address this ultimate case of anisotropic QD 
effects on a LC with strong Q-rj coupling. We per- 
formed complementary optical birefringence and neutron 
diffraction experiments on ?^-dodecyl-l,4-cyanobiphenyl 
(12CB), a rod-like mesogen, spatially confined in colum- 
nar porous silicon (pSi) or porous silica (pSi02) layers. 
These substrates are known to induce anisotropic QD 
O [S]. Both nematic Q and smectic rj OPs are mea- 
sured simultaneously, which allows us to provide unique 
information both about the transition mechanism and 
the hypothetical reduction of the Q-rj coupling. The 
somewhat surprising experimental findings shall be com- 
pared with predictions of state-of-the-art theories of con- 
strained smectic order. 

12CB was purchased from Synthon GmbH and used 
without further purification. This LC compound displays 
the following bulk phase sequence : I for T > 331 K ; SmA 



for 305 K < T < 331 K and crystalline for T < 305 K. 
For neutron scattering experiments, pSi films were pre- 
pared by electrochemical etching of a crystalline, highly 
doped Si wafer [TU]. No special attempt was made to 
remove the native silicon oxide from the inner walls of 
the porous structure. For birefringence measurements, a 
monolithic silica film was prepared in a similar way, fol- 
lowed by thermal oxidation at 800° C for 12 h [U]. Ah 
resulting films exhibit a structure permeated by an array 
of parallel-aligned, non-interconnected channels of 300 
yum length and '^ 10 nm diameter. Transmission electron 
micrographs of the channels |TT] indicate sizeable 1.0 ± 
0.5 nm mean square deviations of their surfaces from an 
ideal cylindrical form. The films were completely filled 
with the LCs by spontaneous imbibition in the isotropic 
phase. Diffraction experiments were carried out on the 
two-axis spectrometer G6.1 at the Leon Brillouin Labo- 
ratory (CEA Saclay). The porous films were macroscopi- 
cally oriented in grazing incidence scattering geometry in 
order to probe momentum transfers parallel to the pores 
axis. The smectic A phase of 12CB is characterized by a 
single observable Bragg peak at qq =0.158 A~^ (see inset 
of Figjll;) . Its integrated intensity is proportional to the 
square of the smectic order parameter, I{T) ex \r]\'^. In 
agreement with measurements on 8CB ^2 , this peak ap- 
pears only for q vectors aligned with the pSi layer surface 
normal, indicating an homogeneous uniaxial macroscopic 
ordering of the smectic stacking along the long axes of the 
pores [T2] . 

Because of its optical transparency and straight pore 
geometry the pSi02 is particularly suitable for optical po- 
larization studies. Accordingly, the orientational molecu- 
lar ordering inside the pores can be precisely evaluated by 
optical birefringence measurements performed in trans- 
mission geometry (13] . To a good approximation, the 
optical birefringence in the nematic or smectic phase is 
proportional to Q [13]. 

FigJT] a and b display the T-dependency of the optical 
birefringence An(T) and integrated intensity I{T) of the 
neutron diffraction smectic peak measured for 12CB in 
the bulk state and under confinement. All measurements 
have been performed upon slow cooling with a rate of 
0.05° C/min. The neutron data are displayed either nor- 
malized to the value at T|^"^ (bulk), or to their maximum 
value (confined samples). As expected for a first order 
I-SmA transition, the two OPs of bulk 12CB simulta- 
neously jump at TjA, with no detectable pretransitional 
effects in the isotropic phase. The magnitude of the jump 
at the transition An{T) « 0.141 well agrees with previ- 
ous studies [TS]. Confinement obviously induces funda- 
mental modifications of the T-behavior of the two OPs. 
Even for temperatures above Tj'^ there exists a weak 
residual birefringence, characteristic of nematic ordering 
along the pore axis. Interestingly, this paranematic order 
induced by confinement is not accompanied by a similar 
parasmectic order (/(T > Tj'^-') = 0). More important. 
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FIG. 1. (Color online) (a) Birefringence and (b) integrated 
intensity of the smectic neutron diffraction peak of 12CB mea- 
sured upon slow cooling, in the bulk (squares) and confined 
in porous Si02 or porous Si (circles), (c) Smectic correlation 
lengths of bulk 12CB (average isotropic (^iso, solid circles) and 
12CB confined in pSi (^y, open circles). Inset: smectic peak 
shape of confined 12CB at T=326 K (experimental resolution 
in filled gray area, the solid line is the best fit according to 
eqnfl] see text). 



the first-order, discontinuous character of the bulk I ~ A 
transition is replaced by a continuous evolution of both 
the orientational (nematic) and translational (smectic) 
ordering as a function of T. Both OPs indicate a transi- 
tion onset temperature Tj'^ which is about 4 K below 
the bulk transition temperature. At lower temperatures, 
in the confined SmA phase, the absolute magnitude of 
An (after the normalization on porosity P) is roughly 
compatible with only 75% of its magnitude in bulk 12CB. 

Interestingly, in the confined SmA phase, smectic or- 
der remains short-ranged, as indicated by a broad Bragg 
peak in Figjlj;. This observation can be attributed to 
the aforementioned susceptibility of the quasi-long range 
smectic order to QD effects introduced by the pore con- 
finement [H 13] . Accordingly, we managed to fit the shape 
of this diffraction peak by a structure factor characteris- 
tic of short-ranged smectic order jl6| while using a partial 



powder averaging procedure, described elsewhere [2]: 
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In Eq. [Tl the first term describes thermally-induced 
smectic fluctuations of dynamical character, while the 
second one accounts for possible static fluctuations in- 
duced by any kind of pore wall irregularities, and thus 
by a QD field [2]- In absence of specific data for 12CB 
in the literature, values for c and for the ratio -A- were 
assumed similar to those of 8CB |2]. We find that the 
broadening of the Bragg peak is entirely attributable to 
the static disorder term at all Ts investigated. An equally 
good fit could be obtained by setting Atherm = 0. This 
differs from 8CB in pSi, which exhibits pre-transitional 
thermal fluctuations reminiscent of the continuous char- 
acter of the N-SmA transition in the bulk ^ [3]. On 
the other hand, it is consistent with the fact that we 
do not observe any pre-transitional smectic diffuse scat- 
tering above 326 K. The temperature dependence of the 
smectic correlation length ^y extracted from the fitting 
procedure is shown in FigjTfc. For confined 12CB, fy rises 
continuously to moderate values upon decreasing T. This 
is in striking contrast with bulk 12CB, where we observe 
a jump of ^ at T^^"^ to resolution limited values. This 
observation provides another strong evidence that the I- 
SmA transition of 12CB in anisotropic pSi actually loses 
its first order character unlike lOCB in isotropic random 
confinement, where a discontinuous variation of ^y was 
retained [3 [9] . It is noteworthy that in case of a strong 
Q-r] coupling, the criticality of the I-SmA phase transi- 
tion might result from the coupling between the nematic 
director and the unidirectional porous geometry, as in the 
case of the I-N transition of 8CB [13] , but also from the 
anisotropic QD effects. These two aspects are discussed 
in what follows. First we examine if this overall new 
phenomenology could be rationalized in the frame of a 
Landau-de Gennes approach. The KKLZ model [HI [TB] 
applied to nCB (n>10) accounts for the Q-rj coupling, 
and also for the anisotropic coupling of the porous geom- 
etry with Q ("field" effect) [If^. The reduced free energy 
then writes : 



=t„Q2 _ 2Q3 + q4 _ ^^Q ^ A<*^g2 
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where i„ and ts are reduced temperatures, cr„ and 
as represent the coupling terms of resp. nematic and 
smectic orderings with the unidirectional field imposed 
by the porous solid, and D is the coupling strength be- 
tween nematic and smectic OPs. While i„ and tg only 
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FIG. 2. (Color online) Equilibrium solutions of the KKLZ 
model in the strong NS coupling limit (_D=2000). All curves 
are normalized at low T, and so that Q{Tj-n) =1 for a hy- 
pothetic system with D = 0. Dashed lines : nematic OP (Q), 
solid lines : smectic OP (r;). Main figure: solutions for a„ = 
0, 1, 5, 15 and 50 and as — 0. Inset: solution for a„ = and 
as = 0, 0.002 and 0.005. 



depend on intrinsic properties of the LC, ct„ and as de- 
pend both on the LC and on the topology of the confining 
matrix [171 llSj . The reduced free energy g is normal- 
ized with respect to nematic ordering, which implies a 
dimensionless prefactor fl weighting the smectic contri- 
bution, il as well as e are material dependent. We select 
il ~ 79000 to map experimental data of the family of 
cyanobiphenyls according to refs. [El [18] . In the KKLZ 
model, disorder effects are simply taken into account by a 
renormalization (lowering) of the transition temperatures 
(tn ->■ ti'" = tn + Ati") [Bj. Fig|2]shows the equihbrium 
values of both nematic Q and smectic 77 OPs as computed 
from the KKLZ model with a Q-rj coupling parameter D 
~ 2000 (strong coupling limit), without disordering terms 
(Ai*^ = 0). Clearly, both (t„ and as field parameters 
induce an increase of the transition temperature. The 
opposite shift observed in the experiment can be quali- 
tatively accounted for by a non-zero disorder term Af!^^ 
in the free energy expansion (Eq. [21 without changing 
the shape of the curves. Above the computed transition 
temperature, we see that the nematic field parameter a^ 
promotes paranematic ordering. This appears as a gen- 
eral feature of unidirectional confinement of LCs [TJl [H] . 
On the other hand, this orientational field contribution 
does not induce similar paiasmectic order, even in this 
strong coupling regime. These two predictions qualita- 
tively agree with experimental observations above Tj°^ , 
provided the field parameter an is rather weak (between 
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FIG. 3. (Color online) Integrated intensity of the smectic 
peak to the power i and birefringence to the power 4 with 
P = 0.17. 



1 and 1.5). These values are in agreement with the ampli- 
tude of field effects attributed to pSi through the analysis 
of the J - TV transition of confined 7CB and 8CB 1^. 
However, the model does not catch the essential charac- 
teristics of the transition itself. Even for very strong cr„ 
values, it rather predicts a remaining first order character 
of the smectic ordering occuring from a paranematic or 
nematic high T-phase. Taking into account an additional 
smectic field term {us ^ 0, see inset of Fig. [2]) the system 
is driven accross a critical threshold, separating continous 
from discontinuous smectic ordering, but it also displays 
a strong parasmectic order at high T, in obvious contrast 
to the experimental observation (see Figlll). 

These discrepancies suggest that random field effects 
play here a dominant role, which is supported by the 
observed downshift of T/^ and the typical Lorentzian 
squared line shape of the diffraction pattern. As expected 
from Garland and lannachione arguments [1] , the transi- 
tion here has crossed a tricritical point under anisotropic 
QD, in a same way as if the Q-rj coupling would be 
turned off f^. The actual continuous character of the 
transition can be checked by inspecting the scalings of 
nematic and smectic OPs as a function of T. As shown 
in FigjSJ both OPs convincingly scale with the power law 
0, ry ex (Tc - TY with T^ « 326.4 K. However, the corre- 
sponding best fit value of /3 = 0.17±0.01 is significantly 
lower than the one expected in the frame of the undis- 
turbed 3D XY model, and which was considered as the 
limit observed for nanoconfined 8CB under QD [H [51 [S] . 
It is rather comparable to temperature scalings observed 
for 8CB in aerogels or pSi [2] ■ The observation of a low 
/? value for a continuous transition is in fact not unusual, 
and is observed in the case of dilute antiferromagnets or 
disordered ferroelectrics [111 HHj : which are prototypical 



realizations of the three dimensional random field Ising 
model (3DRFIM). 

In conclusion, we report combined measurements of 
orientational and translational OPs of a LC confined 
in columnar pSi, characterized in the bulk by a strong 
first order nucleation-growth mechanism. In the high- 
T phase, the confinement induces pa,Ta,nematic but no 
parasmectic order. This fact can be interpreted as re- 
sulting from the unidirectional geometry of pSi, acting 
as an external field coupling to the nematic OP, as in 
the case of the I — N transition of other cyanobiphenyls. 
In the case of 12CB, such field effects are not expected 
to be strong enough to modify the first order character 
of the I — A transition. At Tj°^^ however, no disconti- 
nuity of the OPs nor of the smectic correlation lengths 
is observed. Both Q and 77 grow according to the same 
power law. These observations rather indicate a critical- 
like behavior at low temperatures, which cannot be sim- 
ply accounted for by available phenomenological consid- 
erations. Together with the downshift of T/^ and the 
typical Lorentzian squared line shape of the diffraction 
pattern, they strongly support the prevalence of QD ef- 
fects. Although up to now no available theoretical frame 
can account for the observed transition behaviour, our 
observation of a unique scaling law for both nematic and 
smectic OPs emphasizes that the phenomenology of LCs 
under anisotropic QD cannot be reduced to a simple tun- 
ing of the de Gennes' Q-r] coupling by the strength of the 
random field. 
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